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Harnessing spins as carriers for information has 
emerged as an elegant extension to the transport 
of electrical charges 1 . The coherence of such spin 
transport in spintronic circuits is determined by 
the lifetime of spin excitations and by spin diffu- 
sion. Fermionic quantum gases are a unique sys- 
tem to study the fundamentals of spin transport 
from first principles since interactions can be pre- 
cisely tailored and the dynamics is on time scales 
which are directly observable 2 n . In particular at 
unitarity, spin transport is dictated by diffusion 
and is expected to reach a universal, quantum- 
limited diffusivity on the order of h/m. Here, 
we study the non-equilibrium dynamics of a two- 
dimensional Fermi gas following a quench into 
a metastable, transversely polarized spin state. 
Using the spin-echo technique 12 , we measure the 
yet lowest transverse spin diffusion constant ! 13 1 14 1 
of 0.25(3) h/m. For weak interactions, we observe 
a coherent collective transverse spin-wave mode 
that exhibits mode softening when approaching 
the hydrodynamic regime. 

Studying transport in low- dimensional nanostructures 
has a long and rich history because of its non-trivial fea- 
tures and its relevance for electronic devices. The most 
common case, charge transport, has great technological 
implications and it determines the current-voltage char- 
acteristics of a device. With the development of the field 
of spintronicfl, however, also spin transport has moved 
into the focus of the research interest. Spin transport 
has unique properties, setting it aside from charge trans- 
port: Firstly, the transport of spin polarization is not 
protected by mome ntum conservation and is greatly af- 
fected by scattering^D. Therefore, the question arises 
what is the limiting case of the spin transport coefficients 
when interactions reach the maximum value allowed by 
quantum mechanics? Secondly, unlike charge currents 
(which lead to charge separation and the buildup of an 
electrical field, counteracting the current), spin accumu- 
lation does not induce a counteracting force. 

The main mechanism to even out a non-equilibrium 
magnetization M(r,t) = M(r,t)p(r,t) is spin diffusion, 
which is an overall spin-conserving process. Other, spin 
non- conserving, processes are much slower in ultra cold 
Fermi gas experiments and are neglected hereafter. The 
gradient of the non-equilibrium magnetization VM = 
pVM + MVp drives two distinct spin currents^ The 
first term produces a longitudinal spin current parallel to 
M and the second term induces a transverse spin current 
orthogonal to M. These spin currents are proportional 



* VKTi Bz(x) 



FIG. 1. Quench of a two-dimensional Fermi gas in which all 
atoms were initially prepared in the | \) state, (a), A tt/2- 
pulse prepares the Fermi gas polarized in the ^-direction. 
(b) A magnetic field gradient dB z /dx causes the spins to ac- 
quire different phase-angles <f>{x) in the equatorial plane, (c) 
Collisions cant the spins out of the equatorial plane due to the 
identical spin rotation effect. The acquired projection along 
S z together with the motion of the atoms in the harmonic 
trap impedes rephasing of the spins when the magnetic field 
gradient is reversed. The spin states are illustrated in the 
rotating frame. 



to the longitudinal V\\ and transverse spin diffusivity V±, 
respectively. 

In general, the spin diffusivity in D dimensions behaves 
as V = v Imfp/D, where v is the collision velocity and Z m f p 
is the mean free path. The mean-free path is determined 
by the elastic scattering cross section a between atoms. 
For short-range s-wave interactions at unitarity, the cross 
section attains its maximum value allowed by quantum 
mechanics a ~ V^di^ 1 ? where X^b is the de-Broglie 
wavelength of the colliding particles. At high temper- 
atures in the classical gas regime, the spin diffusion con- 
stant hence scales as V oc T D I 2 (up to logarithmic correc- 
tions in two dimensions owing to the energy dependence 
of the cross-sectiorf^) . In the degenerate regime (of any 
dimensionality), the situation is quite different: the de- 
Broglie wavelength X^b is of order 1/kp and hence the 
mean- free path is Z m f p = l/(na) ~ 1/fcp, where kp is the 
Fermi wave vector of the gas and n ~ kp is the density. 
Hence, the spin diffusion constant is given by Planck's 
constant H divided by the mass m of the particle car- 
rying the spin. This quantum limit can also be viewed 



as a result of the uncertainty principle by noticing that 
the mean-free path is limited by the mean interparticle 
spacing^. The simple scaling argument, however, hides 
much of the rich underlying physic s. In par ticular, it can- 
not explain the Leggett-Rice effectP^ 4 ^', the difference 
between longitudinal and transverse spin diffusivitieJ^, 
and the transition to weak interactions where the physics 
changes because the system evolves from hydrodynamic 
to collisionless. So far, lowest spin diffusion constant 
has been measured for longitudinal spin currents to be 
D\\ = 6.3h/m in three-dimensional degenerate Fermi 
gases at unitarity 5 , approximately two orders of mag- 
nitude smaller than in semiconductor nanostructures^. 

Here, we study the coherence properties of a trans- 
versely polarized two-dimensional Fermi gas under the 
influence of magnetic field gradient. By employing the 
spin-echo technique 12 we gain access to the intriguing, 
but little understood, transverse spin dynamics in two- 
dimensional Fermi systems^. The experiment starts 
with a polarised Fermi gas, where all spins are in the 
same state \\). Three resonant rf-pulses are applied to 
rotate the spin state by respective angles tt/2 - it - tt/2 
about a fixed axis in spin space, S x . Consecutive pulses 
are separated by a time r and we refer to the spin evolu- 
tion time as t = It. The initial 7r/2-pulse creates trans- 
verse spin polarisation in a coherent superposition be- 
tween |^) and |t) (see Fig. la). An applied magnetic 
field gradient B' = dB z /dx gives rise to a transverse 
spin wave (see Figure lb) of wave vector Q = 5jtB\ 
where £7 = 152kHz/G is the differential gyromagnetic 
ratio between |t) and \\) at 209 G and t is the evolution 
time (see Methods). On time scales much shorter than 
the trapping period, Q is independent of the interaction 
strength^, and nearby atoms acquire a relative phase- 
angle of A<fi « Q/kp- This lifts the spin-polarisation of 
the Fermi gas and the spins can collide with each other 
as they move in the harmonic potential or diffuse. Triv- 
ial dephasing induced by the magnetic field gradient is 
reversed by the 7r-pulse and the spin state refocus after a 
time r if no decoherence has occured. The final 7r/2-pulse 
maps the spin state from the S y onto S z , which is mea- 
sured by performing a Stern-Gerlach experiment in time 
of flight, and we record (M z ) = (JV t - iY;)/(iV t + N±). 

Even though the interaction potential between atoms 
is not explicitly spin-dependent, an effective spin- 
exchange interaction is mediated by the required anti- 
symmetrisation of the scattering wave function. In bi- 
nary collisions, this leads to the identical spin rotation 
effect 20 : in a collision both spins rotate about the axis 
defined by the sum of their spin orientations. In the 
Fermi degenerate regime, the binary collision picture 
has to be replaced by Landau's quasiparticle description. 
Here, quasiparticle excitations can be considered being 
affected by a "molecular field " resultin g from the effec- 
tive spin-exchange interactio n 1 13 * 14 * 21 -^. As a result, the 
spin wave cants out of the S x — 5^-plane and forms a 
spin spiral which acquires a component along S z (see 
Figure lc). The magnitude of the identical spin rota- 
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FIG. 2. Spin-echo signals in the strongly interacting regime. 
a, Spin echo-signal at \n(kFa>2D) = —0.2. The blue line is a 
fit oc exp[— (2rr) 3 ]. b, Transverse spin diffusion constant T>± 
as a function of interaction strength deduced from the decay 
constant T of the spin-echo signal. Illustration of the different 
spatial variation of Fermi surfaces of a polarized Fermi gas for 
the case of longitudinal (c) and transverse (d) spin diffusion. 



tion effect is proportional to the mean-field interaction 
strength g = —27rh 2 /[m\n(kFCL2D)} of the gas. Here, kp 
is the Fermi wave vector for the initially spin-polarized 
sample, and a2D is the two-dimensional scattering length, 
which is linked to the binding energy of the confinement- 
induced dimer by^i^g = h 2 /ma 2 D . 

In the strongly interacting regime, i.e., —1 < 
lii(kFd2D) < 3, spin transport is dominated by diffu- 
sion and the spin-echo signal attains a characteristic non- 
exponential decay of the formal!] ^ exp| — (2Fr) 3 ]. Fit- 
ting the envelope of the spin-echo signal 24 (see Figure 
2a) with this function we deduce the transverse spin 
diffusion constant from the measured decay rate a d 12 * 25 * 
V± = 12h 2 T 3 /(Sj B f ) 2 . We extract the transverse spin 
diffusion constant for various interaction strengths and 
find V± = 0.25(3)h/m at a shallow minimum around 
ln(fc j pa2i)) = (see Figure 2b). From the arguments 
given at the beginning, we expect the diffusivity to be 
on the order of h/m. Observing a smaller value than for 
longitudinal spin diffusion in three dimensions 5 is prob- 
ably less linked to the dimensionality rather than to the 
phase space available for collisions necessary to drive spin 
diffusion^ i n the case of longitudinal spin currents, a 
gradient of the M z (x) polarization along the x-direction 
can be considered as a spatial variation of the local Fermi 
surfaces kp^(x) and kp,i(x) (see Fig. 2c). Only |t) spins 
in a small region near the Fermi surface can diffuse from 
x to x + dx, invoking the typical T~ 2 scaling for quasipar- 
ticles in the deeply degenerate regime. In contrast, in the 
case of transverse spin currents the Fermi surfaces at dif- 
ferent positions are of the same size but have slightly dif- 
ferent directions of magnetization, (see Fig. 2d). Hence, 
a spin moving along x from anywhere between the Fermi 
surfaces for spin-up and down must scatter to reach lo- 
cal equilibrium, which scales as n^ — n± and provides a 
much larger phase space. The result is that for a degen- 
erate system the transverse diffusivity is smaller than the 
longitudinal and becomes independent of temperature^. 

For a repulsively interacting Fermi gas in three dimen- 
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FIG. 3. Observation of spin wave for the non-interacting 
gas, i.e., ln(/cF«2D) — > +oo: a, We plot the envelope of the 
spin-echo signal (dots) and show the fit function (grey line). 
b, centre of mass separation of |t) and \l) after 11.8 ms of 
time of flight, c and d, momentum distribution of local spin 
polarisation for large and small spin polarisation, respectively. 
The lower two panels in d show the distribution of spin up 
and down, respectively. 



sions the decay of a spin-spiral state into a Stoner itin- 
erant ferromagnet was considered^. In this case, the 
growth rate of the instability towards the formation of 
the ferromagnet order was found to scale ex Q 2 and hence 
the spin-polarisation signal would decay ex exp[— (at) 3 ], 
similar to spin diffusion. 

In the weakly interacting attractive Fermi liquid 
regime 26 , i.e., \n(kF(i2D) > 3, we observe a qualitatively 
different behaviour (see Figure 3a). The envelope of the 
spin-echo signal decays exponentially and is modulated 
by a slow oscillation and both the frequency and the 
damping constant of the oscillation depend on the in- 
teraction strength. We attribute this slow periodic mod- 
ulation to transverse spin waves which are excited by 
the spin-echo sequence. The spin current induced by the 
magnetic field gradient can only be inverted completely 
by the 7r— pulse when applied in phase with the harmonic 
motion. This can be qualitatively (for zero interactions) 
understood in a phase-space picture by considering that 
the magnetic field gradient causes a spatial offset between 
the potentials for the |t) and the \\) spin states and, 
correspondingly, the phase-space trajectories of the two 
spin-components are displacecP. Indeed, we observe that 
the difference of the center-of-mass momenta (k^) — (k±) 
of the two spin densities n^(k) and n^(k) oscillates out- 
of-phase with the contrast of the spin-echo signal (Figure 
3a and 3e). The total density profile n(k) = n^(k)+ni(k) 
is stationary, which indicates a pure spin mode. 

The measured envelope of the spin-echo signal (see 
Fig. 3a) is fitted with an empirical function of the form: 
Aex-p(--yt) (l-B\sm(ujt/2)\ 3 ) + C. Here, \A\,\B\ ~ 
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FIG. 4. Spin wave in the collisionless regime as function of in- 
teraction strength, a, Time evolution of global spin polariza- 
tion for different interaction strengths, b, Frequency of oscil- 
lating spin-polarization. For the data point at ln(/cF«2D = 3) 
we were not able to fit a frequency and set the value to zero. 
c Damping coefficient 7 of the spin-polarization. The dashed 
line indicates the elastic scattering rate and the solid line is a 
fit to the measured decay. 



(9(1) and \C\ <C 1 are amplitudes and a global offset 
of the spin echo signal, respectively. We find that this 
periodic function of frequency uj is approximating the 
non-sinusoidal signal best. The exponent is fixed to the 
value 3 for all interactions, and we do not see a large in- 
fluence on the deduced parameters. Both the frequency 
uj and the exponential damping 7 depend on the interac- 
tion strength of the gas. 



For zero interaction strength the mode is at half the 
trap frequency and we observe a softening of the spin 
mode when approaching the hydrodynamic regime (Fig- 
ure 4a and b). Such a behaviour is reminiscent of 
the collisionless-to-hydrodynamic transition of the spin- 
dipole oscillation^^, whereas spin-symmetric modes 
have non-zero frequencies even in the hydrodynamic 
regime 28 . The decay rate 7 of the oscillation (see Figure 
4c) is found smallest for the non-interacting gas, where 
it is two orders of magnitude smaller than the dephasing 
rate l/rR amsey in the inhomogeneous magnetic field (see 
Methods). As the interaction strength is increased, the 
rate 7 is growing (see Fig. 4c). We fit this growth with 
a power-law behaviour 7 = a/\n(kFQ>2D)^ + 7o and we 
determine the exponent to be j3 = 1.3 ± 0.15. Observ- 
ing a scaling with an exponent close to unity signals that 
the loss of global coherence is dominated by the identi- 
cal spin-rotation effect, which scales proportional to the 
mean-field interaction strength g. In contrast, the rate 
of elastic collisions scales as ex l/ln^F^Z)) 2 , i-e., (3 = 2 
(see dashed line in Fig. 4c). This suggest the following 
decoherence mechanism in a spin- 1/2 Fermi gas: after 
building up in the magnetic field gradient, the spin wave 
decays when neighbouring atoms have acquired a suffi- 
ciently large phase-angle Ac/) between their spins such 
that in a spin-conserving collision both spins rotate out 
of the S x — 5^-plane. Initially, this is a coherent pro- 
cess, which subsequently dephases due to the coupling of 
the spin-rotation to the momentum of the atoms 29 . This 
picture contrasts to the previously conjectured mecha- 
nism of the decoherence of a sp in mixture resulting from 
a magnetic field curvatur e! 3 1 30 1 , which has been probed, 
for example, using the Ramsey technique^. Our spin- 
echo measurements show that the dephasing time mea- 
sured by the Ramsey technique underestimates than the 
actual spin decoherence time in the weakly interacting 
regime by orders of magnitude, affecting the question of 
the equilibration of the Fermi gas. 
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I. METHODS 

A. Preparation of the two-dimensional gases 

We prepare a quantum degenerate Fermi gas of 40 K 
atoms in a one-dimensional optical lattice of wavelength 



A = 1064 nm populating a stack of approximately 40 in- 
dividual two-dimensional quantum gase J 31 | 32 4 The radial 
confinement of the two-dimensional gases is harmonic 
with a trap frequency of uj r = 2tt x 127 Hz and the ax- 
ial trap frequency is u z = 2n x 75 kHz. Using a radio- 
frequency cleaning pulse, we ensure preparation of a spin- 
polarized gas in the \F = 9/2, mp = —9/2) = |^) ground 
state with 1.5 x 10 5 atoms, corresponding to the Fermi 
energy of the spin-polarised gas Ep = 12.3(2.0) kHz at 
a temperature of ksT/Ep = 0.24(3). The coupling 
strength of a two-dimensional Fermi gas is given by the 
parameter g = — 27rh 2 /[m\n(kFd2D)r^- Here, kp is the 
Fermi wave vector for the initially spin-polarized sample, 
and d2D is the two-dimensional scattering length which is 
linked to the binding energy of the confinement-induced 
dimer by Eb = h 2 /ma 2 D . 

Our experiments on the spin dynamics start by ramp- 
ing the magnetic field from B = 209.15 G to the desired 
value B near the Feshbach resonance in 30 ms and wait 
for 200 ms. Subsequently we apply the spin-echo pulse 
sequence (see main text) of three radio-frequency pulses 
with frequencies of ~ 50 MHz. The pulses have a square 
envelope and duration of t^/2 ~ 46 /is for a 7r/2-pulse and 
2 x t^/2 f° r a 7r-pulse. The pulses are separated by the 
same time r which gives the total spin evolutions time 
t = 2r. 



Ramsey measurements: Dephasing in the 
magnetic field gradient 



We assess the time scale of simple dephasing of the 
spin in the magnetic field gradient by performing Ramsey 
spectroscopy on the non-interacting gas using the follow- 
ing sequence: the initial 7r/2-pulse is followed by a sec- 
ond 7r/2-pulse with a variable time delay r. We observe a 
Ramsey coherence time of TR arnsey = (600±30)/is. Since 
the Ramsey time is more than an order of magnitude 
shorter than the inverse trap frequency, we neglect the 
motional contribution to the dephasing. Instead, we re- 
late the Ramsey time to the dephasing time of the (quasi 
stationary) spins across the whole cloud, and we find from 
this the magnetic field gradient. We fit the data with 
P{t) = 8NJ ^ W Where, = 2irAt = d^B'R F t. A is the 
Zeeman shift due to the gradient at the edge of the cloud 
compared to the center. The fit function fits <j) = bt, i.e., 
b = 2ttA. The fit results gives b = 8.1(1) kHz and we 
know dj = 152kHz/G at 209.15 G and R F = 17.7/im. 
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